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Abstract 

In this paper we introduce an algebraic form of the dispersion relation for a non 
ideal inhomogeneous dusty plasma in order to improve drastically the calculation of 
the drift instability growth rate. This method makes use of the multipole approx- 
imation of the Z dispersion function, previously published, and valid for the entire 
range. A careful analysis of the solutions spectra of this kind of polynomial equation 
permits us to calculate easily the growth rate of the drift instability for the ion-dust 
and dust acoustic mode. The value of the parallel to magnetic field wavelength for 
which the instability reaches the maximal value is carefully localized and discussed. 
The unstable dust-ion and dust acoustic mode are discriminated and analyzed in 
function of the density gradient, Te/Ti - ratio, and dust grain radius. 
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1 Introduction 

Plasma inhomogeneities across the magnetic field in the presence of finite - size charged 
grains causes a wide class of instabilities of an inhomogeneous dusty plasma called gradient 
instabilities. Such instabilities can be studied in the approximation on magnetic field 
where we have parallel straight field lines in order to simplify our treatment. We look 
for instabilities in the very low frequency regime where a new spectrum instabilities and 
waves appear, induced by the dust collective dynamics: Dust - Acoustic - Waves (DAWs), 
Dust - Ion - Acoustic - Waves (DIAWs), etc. The frequency of DAWs are around 10 Hz 
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as determined in the laboratory and lower in astrophysical plasmas [1,2]. In the case that 
grains are in the micron range we expect a non - ideal behavior due to the fact that the 
particulate are highly charged and intermolecular forces could play certainly an important 
role. In order to discuss this problem we compare the ideal properties with the simple 
hard - core model and in a next work we will use a better model by means of of the square 
- well model and the Pade rational approximant to the equation of state [3] for hard - 
sphere gas, that in our knowledge is more realistic as the simple application of the Van 
der Waals equation of state [4]. In this paper we show an analysis of the electrostatic 
waves and instabilities growth rates in a weakly non - ideal magnetized dusty plasma with 
density and temperature gradients, ignoring charge fluctuation. As introduced before, the 
non - ideal behavior is characterized by the hardcore model defined by 

p = nT(l + b Q n), 

or in similar manner by the square - well model given by the Ree and Hoover expression 
[5]. 



2 Theoretical Model 

In this paper we introduce a new numerical treatment in combination with a more realistic 
formulation of the equation of state to simulate weak non ideal effects in order to analyze 
inhomogeneous Vlasov - Dusty Plasma systems where a linearized dispersion relation 
is obtained. Due to the lower frequency range (u,k z VT <C w c ), enough energy can be 
transferred from the particle to the wave and instabilities can be generated. In order to 
get an adequate linear dispersion relation with a magnetic field given by B = B^k for 
Maxwellian multi-species plasmas (electron, ion and dust), we introduce our well known 
and very accurate multipolar approximation [6] for the Z dispersion function. 
In the presence of a magnetic field we have the distribution function of the species a, 
solution for the kinetic equation 

dfa _ , df oa , . 

-TT = V0-— — 1 

at m a av 
in the time dependent following form [7, 8] 

/(r, y,t) = ^ f exp [iu{t - t')\ V0'(r(t')) • £f±rdt (2) 

where a — e, i, d. Now, the dispersion relation in terms of the dielectric susceptibilities, 
in the low frequency approximation (u>, k z Vx u> c ) is 

l + $> oa = (3) 
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where, 
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Further, in order to simplify our expressions, we use: 
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Now, using the following identity for the dispersion function Z 



(5) 



Z' = -2[l + £ a Z(£ a )], 



we obtain after several cumbersome algebraic manipulations the dielectric susceptibility 
in the form 
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In order to put our dispersion relation in a dimensionless form, we introduce following 
suitable definitions: 
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Now, using those results and assuming that uj <^ uj oi <^ uj 0( i we can write down Eq.(3) as 

1 + XOe + XOi + X0d = (7) 
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In the non ideal case (dust) we introduce a relation that in principle express the non ideal 
behavior of the system in terms of the pressure in the form 



p = n° d T d (l + b d n° d ) 



(8) 



given by the hard-core model. This model is taken for simplicity. A better model, as 
mentioned before, will be introduced in a future work. Now, following definitions are also 
useful 

1 Vp d 1 V< 1 VT d 
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Those relations are very convenient by writing the full dispersion relation[4]. In fact we 
have 

1 _ 1 + 2b d n 0d 1 | 1 ^ 
L p 1 + b d n od L n L T 

for the non-ideal case. For the ideal one, we use the well known relation p Q j = n jTj, and 
in a similar way we get 

^ = ^ + ^ (Hi 

Lpj Lrij L Tj 

where j = i, e. Two special cases can be worked out: 

A) Density gradient equal to zero Vnoj = 0, that means, Lpj = L^- 



B) Temperature gradient equal to zero VI} = 0, that means, Lpj = Lrij. 

Further we can introduce following relations in order to express dielectric susceptibilities 

in a suitable forms 
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Using those relations we arrive to the dispersion relation for the case B where we get: 
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where A„ d = [(1 + 2b d n 0d )/ (1 + b d n od )]A p and A n = A 



Pi ■ 
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In a similar way, it is possible to include the terms for case A, where we shall have 

A nj =A Pr (17) 

Introducing now the multipolar approximation to Z we can get a polynomial expression 
in the well known form [9] 

^aiW/^bjW = (18) 

* 3 

where coefficients Oj and bi are functions of the system parameters. Such an expression is 
easy to solve and with high accuracy to find roots of the numerator. An analysis of these 
solutions spectra permit us to give the imaginary parts 7 = Im{VL) in function of 1/K y , 
which represent the growth rate instabilities. 



3 Results and Conclusions 

The quasi-neutrality equation for dusty plasmas can be approached by a simplified one 
due to the high state of charge of the dust grains 

n oi = Z D n od + n oe ~ Z D n od (19) 

and the electron susceptibility can be neglected in the dispersion relation. The range 
of the main parameters in the study of the low frequency oscillation of dust grains is 
established by the approximations that conduced to the simplified dispersion relation 

n, K z u d < n cd (20) 

Unstable dust oscillations (Im(Q) > 0) are found for Vt cd ~ 10 _1 , K z U d ~ 1CT 2 . At 
the present time, we only give the results for the density gradient case {i.e. d/dT = 
0). For slightly inhomogeneous plasmas with normalized density gradient length A n = 
n oa / (Az?iVn oa ) ~ 10 2 , the shape of the dust instability (Jm(fi) max ) curve as function of the 
perpendicular to magnetic field wavelength (1/K y ) is similar to that for ions, previously 
studied [8]. 

The maximum value of the instability increases and narrows with the state of charge of 
the dust Zd but decreases and get wider with the mass. For typical laboratory light dusty 
plasmas (m d ~ 10 4 m p , Z D ~ 10 3 ) the instability of dust acoustic or electrostatic waves is 
narrower and smaller than that for ions In figure 1 the peak of the left corresponds to the 
typical shape of instability of slightly inhomogeneous plasmas, while the right region of 
instability appears for density gradient lengths of the order of a hundred of Debye lengths 
(A n = A < 10 2 ). For higher density gradients (A = 5 x 10 1 ), this new instability region 
is wider and so high as the typical one. For even higher density gradients (A = 10 1 ), 
figure 2 shows that the new right region gives a higher instability. This figure also shows 
the effect of the non ideality of the plasma, necessary condition for the exhibition of 
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Figure 1: Normalized maximum growth rate as a function of normalized perpendicular 
wavelength for slightly inhomogeneous plasma (A = 10 2 ) and for a relatively inhomoge- 
neous one (A = 5 x 10 1 ). 
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Figure 2: Normalized maximum growth rate as a function of normalized perpendicular 
wavelength for ideal and non ideal plasmas. 

dust acoustic waves. For typical laboratory dust radius the b Q d parameter of the hard 
core potential equation of state, is of the order of 10~ 14 m 3 . And for typical values of 
ion density of 10 16 m -3 (and corresponding n d, by quasi-neutrality relation), it appears 
a new intermediate instability region which can reach a maximum for denser plasmas 
(n i = 10 17 m~ 3 ) or larger dust particles (b d > 10~ 13 m 3 ). This maximum is limited for 
the relation for dust collective behavior 

r d < \ Di (21) 
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